Progress in Electromagnetic Research, M (PIER M, 14, (2010), 193-206. 



O 

(N 

C 

a 
m 

(N 



Oh 



(N 
> 
in 

(N 



i 



Electromagnetic wave scattering by many small 
particles and creating materials with a desired 

permeability 

A.G. Ramm 

Department of Mathematics 
Kansas State University, Manhattan, KS 66506-2602, USA 

rammOmath . ksu . edu 
Abstract 

Scattering of electromagnetic (EM) waves by many small imedance 
particles (bodies), embedded in a homogeneous medium, is studied. 
Physical properties of the particles are described by their boundary 
impedances. The limiting equation is obtained for the effective EM 
field in the limiting medium, in the limit a — > 0, where a is the charac- 
teristic size of a particle and the number M(a) of the particles tends 
to infinity at a suitable rate. The proposed theory allows one to cre- 
ate a medium with a desirable spatially inhomogeneous permeability. 
The main new physical result is the explicit analytical formula for the 
permeability fi(x) of the limiting medium. While the initial medium 
has a constant permeability /io, the limiting medium, obtained as a 
result of embedding many small particles with prescribed boundary 
impedances, has a non-homogeneous permeability which is expressed 
analytically in terms of the density of the distribution of the small 
particles and their boundary impedances. Therefore, a new physical 
phenomenon is predicted theoretically, namely, appearance of a spa- 
tially inhomogeneous permeability as a result of embedding of many 
small particles whose physical properties are described by their bound- 
ary impedances. 

PACS: 02.30.Rz; 02.30.Mv; 41.20.Jb 
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Key words: electromagnetic waves; wave scattering by many small bodies; 
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1 Introduction 



In this paper we outline a theory of electromagnetic (EM) wave scattering 
by many small impedance particles (bodies) embedded in a homogeneous 
medium which is described by the constant permittivity eo > 0, permeabil- 
ity fio > and, possibly, constant conductivity o"o > 0. The small particles 
are embedded in a finite domain 0. The medium, created by the embed- 
ding of the small particles, has new physical properties. In particular, it 
has a spatially inhomogeneous magnetic permeability /u(x), which can be 
controlled by the choice of the boundary impedances of the embedded small 
particles and their distribution density. This is a new physical effect, as far 
as the author knows. An analytic formula for the permeability of the new 
medium is derived: 

where 

^j(x) = l-| h(x)N(x). 

Here uj is the frequency of the EM field, /io is the constant permeability 
parameter of the original medium, h(x) is a function describing boundary 
impedances of the small embedded particles, and N(x) > is a function 
describing the distribution of these particles. We assume that in any sub- 
domain A, the number Af(A) of the embedded particles D m is given by the 
formula: 

A/"(A) = [ N(x)dx[l + o(l)], a^O, 

a 2 K J a 

where N(x) > is a continuous function, vanishing outside of the finite 
domain f2 in which small particles (bodies) D m are distributed, k G (0, 1) is 
a number one can choose at will, and the boundary impedances of the small 
particles are defined by the formula 

_ h{x m ) 
a K 

where point inside m— th particle D m , Re h(x) > 0, and h(x) is a 

continuous function vanishing outside f2. The impedance boundary condition 
on the surface S m of the m— th particle D m is E l = (^ rn [H t , N], where E l 
(H 1 ) is the tangential component of E (H) on S m , and N is the unit normal 
to S m , pointing out of D m . 

Since one can choose the functions N(x) and h(x), one can create a 
desired magnetic permeability in O. This is a novel idea, to the author's 
knowledge, see also [TT] . 
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We also derive an analytic formula for the refraction coefficient of the 
medium in created by the embedding of many small particles. An equation 
for the EM field in the limiting medium is derived. This medium is created 
when the size a of small particles tends to zero while the total number 
M = M(a) of the particles tends to infinity at a suitable rate. 

The refraction coefficient in the limiting medium is spatially inhomoge- 
neous. 

Our theory may be viewed as a " homogenization theory", but it dif- 
fers from the usual homogenization theory (see, e.g., pQ, [5], and references 
therein) in several respects: we do not assume any periodic structure in the 
distribution of small bodies, our operators are non-selfadjoint, the spectrum 
of these operators is not discrete, etc. Our ideas, methods, and techniques 
are quite different from the usual methods. These ideas are similar to the 
ideas developed in papers 0[8], where scalar wave scattering by small bodies 
was studied, and in the papers [9]. [10). However, the scattering of EM waves 
brought new technical difficulties which are resolved in this paper. The dif- 
ficulties come from the vectorial nature of the boundary conditions. Our 
approach is valid for small particles of arbitrary shapes, but for simplicity 
we assume that the small bodies are balls of radius a. 

We give a new numerical method for solving many-body wave-scattering 
problems for small scatterers. 

In Section 2 an outline of our theory is given and the basic results are 
formulated and explained. In Section 3 the Conclusions are formulated. In 
Section 4 short proofs of two lemmas are given. In Section 5, Appendix, 
more difficult and lengthy proofs are given, and a numerical method for 
solving many-body wave scattering problem in the case of small scatterers 
is presented. 

2 EM wave scattering by many small particles 

We assume that many small bodies D m , 1 < m < M, are embedded in a 
homogeneous medium with constant parameters eo, fJ-o- Let k 2 = u 2 eofio, 
where uj is the frequency. Our arguments remain valid if one assumes that 
the medium has a constant conductivity do > 0. In this case eo is replaced 
by eo + i^j- Denote by [E,H] = E x H the cross product of two vectors, 
and by (E, H) = E ■ H the dot product of two vectors. 

Electromagnetic (EM) wave scattering problem consists of finding vec- 
tors E and H satisfying the Maxwell equations: 

V x E = icjfj, H, V x H = -iu;e Q E in D := R 3 \ U^ =1 D m , (1) 
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the impedance boundary conditions: 

[TV, [E, N}} = ( m [H, N] onS m , l<m< M, (2) 

and the radiation conditions: 

E = E + v E , H = H + v H , (3) 

where Cm is the impedance, N is the unit normal to S m pointing out of 
D m , Eq,Hq are the incident fields satisfying equations (pQ) in all of M 3 . One 
often assumes that the incident wave is a plane wave, i.e., Eq = £e lkax , £ 
is a constant vector, a S S 2 is a unit vector, S 2 is the unit sphere in R 3 , 
a ■ £ = 0, ve and vh satisfy the radiation condition: r(|^ — ikv ) = o(l) as 
r := \x\ — > oo. 

By impedance Cm we assume in this paper either a constant, Re Cm > 0, 
or a matrix function 2x2 acting on the tangential to S m vector fields, such 
that 

Re(C m M^)>0 V£*eT m , (4) 

where T m is the set of all tangential to »S m continuous vector fields such 
that DivE 1 * = 0, where Div is the surface divergence, and E is the tan- 
gential component of E. Smallness of D m means that ka <C 1, where 
a = 0.5 maxi< m <M diamD m . By the tangential to S m component E t of 
a vector field E the following is understood in this paper: 

E l = E - N(E, N) = [N, [E, N}}. (5) 

This definition differs from the one used often in the literature, namely, from 
the definition E l = [N, E]. Our definition ([5]) corresponds to the geometrical 
meaning of the tangential component of E and, therefore, should be used. 
The impedance boundary condition is written usually as 

E^CiH^N], 

where the impedance <C is a number. If one uses definition (|5|), then this 
condition reduces to (|5J), because [[TV, [H,N]],N] = [H,N\. 

Lemma 1. Problem h> as a t most one solution. 

Lemma [1] is proved in Section 2. 
Let us note that problem (JT|)- (SJ) is equivalent to the problems ([6]), (JJJ)) ([3]), 
@, where 

V x V x E = k 2 E in D, H = V X E , (6) 
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[N, [E, N}} = -^-[V x E, N] on S m , l<m< M. (7) 

Thus, we have reduced our problem to finding one vector E(x). If E(x) is 
found, then H = 

Let us look for E of the form 

E = E +Y^V x / g(x, t)a m {t)dt, g(x, y) = — (8) 

, is 47tx — y\ 

m=l J b ™ 

where t 6 S m and dt is an element of the area of S m , a m (t) € T m . This E 1 
for any continuous a m {t) solves equation ([6]) in D because Eq solves ([6]) and 

VxVxVx / g(x,t)a m (t)dt = VV • V x / g(x,t)a m {t)dt 

- V 2 V x / g(x,t)a m (t)dt 

J Sm 

= k 2 V x / g(x,t)a m (t)dt, x £ D. 

J S m 

(9) 

Here we have used the known identity divcurlE = 0, valid for any smooth 
vector field E, and the known formula 

- V 2 g(x,y) = k 2 g{x,y) + 5{x-y). (10) 

The integral f s g(x,t)a m (t)dt satisfies the radiation condition. Thus, for- 
mula ([5D solves problem ([7|), ©, 0D, if cr m (t) are chosen so that bound- 
ary conditions ([7]) are satisfied. 

Define the effective field E e (x) = E™(x) = Ei m \x,a), acting on the 
m— th body D m : 

E e (x) := E{x) - V x f g(x,t)a m (t)dt:=Ei m \x), (11) 

where we assume that x is in a neigborhood of S m , but E e {x) is defined 
for all x € M 3 . Let x m £ D m be a point inside D m , and d = d(a) be the 
distance between two neighboring small bodies. We assume that 

lim-^ = 0, limd(a) = 0. (12) 
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We will prove later that E e (x,a) tends to a limit E e (x) as a — > 0, and 
E e (x) is a twice continuously differential) le function. To derive an integral 
equation for a m = a m (t) , substitute 



E = E e + V x I g(x,t)a m (t)dt 



into impedance boundary condition ([7j), use the known formula (see, e.g., 

my- 



[N,V x j g(x,t)a m (t)dt] T = j [N s ,[S7 x g(x,t)\ x=s ,a m (t)]}dt ± 



2 ' 

(13) 

where the ± signs denote the limiting values of the left-hand side of as 
x — > s from D (D m ), and get the following equation: 

a m (t) = A m a m + f m , 1 < m < M. (14) 

Here A m is a linear Fredholm-type integral operator, defined by formula 
A m = — 2[N S , B m arn\, where the operator B m is defined by formula ([20]) . 
and f m is a continuously differentiable vector function, defined by formula 

(USD. 

Let us find formulas for A m and f m . Equation (|14p is derived in Ap- 
pendix and there the formulas for f m and A m are obtained. 
One has: 

f m = 2[f e (s),N s ], f e (s) := [N s ,[E e (s),N s ]]-^[V xE e ,N s }. (15) 
Boundary condition (J7J) and formula (|13[) yield 

fe(s) + ~[<r m (s),N a ] + [ j [N s , [V s g(s,t),a m (t)]]dt,N s 

Cm 



(16) 



[V x V x / g(x,t)a m (t)dt, N s ] \ x ^ s = 0. 
^Mo Js m 

Using the known formula V x Vx = graddiv — V 2 , the relation 
V^V^ • / g(x,t)a m (t)dt = V X / (-V t g(x,t),a m (t))dt 

J Sm " Sm 

= V X g(x,t)Biva m (t)dt = 0, 

J Sm 



(17) 
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where Div is the surface divergence, and the formula 

-Vl I g(x,t)a m (t)dt = k 2 I g(x,t)a m (t)dt, x G D, (18) 

where equation pup was used, one gets from (|16p the following equation 

- [N„ a m (s)] + 2f e (s) + 2Ba m = 0. (19) 

Here 



Ba m := [ / [N s , [V s g(s,t),a m (t)]}dt,N s }+( m iujeo[ / g(s,t)a m (t)dt,N s ]. 

(20) 

Take cross product of N s with the left-hand side of (|19|) and use the formulas 
N s ■ a m (s) = 0, f m := f m (s) := 2[f e (s),N s ], and 

[N s ,[N s ,a m {s)}} = -a m {s), (21) 

to get from (fT9|) equation (JT3J): 

(Tm(s) = 2[f e {s),N s ] - 2[N s ,Ba m ] := A m a m + f m , (22) 

where 

A m cr m = ~2[N S , Ba m ]. 

The operator A m is linear and compact in the space C(S m ), so that equation 
(j22j) is of Fredholm type. Therefore, equation (j22|) is solvable for any f m € 
T m if the homogeneous version of (I22p has only the trivial solution a m = 0. 
In this case the solution a m to equation (|22p is of the order of the right-hand 
side f m , that is, 0(a~ K ) as a —> 0, see formula (]15p . Moreover, it follows 
from equation (|22p that the main term of the asymptotics of a m as a — > 
does not depend on s € # m . 

Lemma 2. Assume that a m G T m , cr m G C(S m ), and a m (s) = A m a m . Then 
a m = 0. 

Lemma [2] is proved in Section 2. 

Let us assume that in any subdomain A, the number A/"(A) of the em- 
bedded bodies D m is given by the formula: 

Af(A) = ^—[ N(x)dx[l + o(l)l, a^0, (23) 
a 2 K J A 
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where N(x) > is a continuous function, vanishing outside of a finite domain 
in which small bodies D m are distributed, n £ (0, 1) is a number one can 
choose at will. We also assume that 

_ h{x m ) 
a K 

where Re h(x) > 0, and h(x) is a continuous function vanishing outside 0. 
Let us write (|8|) as 

M M . 

E(x) = E (x)+'^2[V x g(x,x m ),Q m \+'^2 Vx / (g(x,t)-g(x,x m ))a m (t)dt, 

m=l m=l JSm 

(25) 

where 

Q m := I u m (t)dt. (26) 

The central physical idea of the theory, developed in this paper, is simple: 
one can neglect the second sum in (|25p compared with the first sum. 

Since a m = 0(a _K ), one has Q m = 0(a 2 ~ K ). We want to prove that the 
second sum in (|25p is negligible compared with the first one. This proof is 
based on several estimates. 

We assume in these estimates that a — > 0, | — >• 0, and |cc — x m | > d. 
Under these assumptions one has 

ji ■= \[V x g(x,x m ),Q m )\ < O (maxjl ^V) °(« 2 ^) 5 ( 27 ) 

32 ■= |V x jf (s(x,t) -5(x,x m ))<r m (t)dt| < aO fmax | , ^ | ^ 0(a 2 ~ K ), 

(28) 

and 

j = O (max{^,/ca}) -»• 0, ^ = o(l), a -> 0. (29) 

These estimates show that one may neglect the second sum in (|25p . and 
write 

M 

E e {x) = E (x) + [Vx9(x, x m ),Q m ] (30) 

m=l 

with an error that tends to zero under our assumptions as a — > 0, and when 
\x — Xj\ ~ a then the term with m = j in the sum (|30p should be dropped 
according to the definition of the effective field. We will show that the limit 
of the effective field, as a — » does exist and solves equation (|39l) . see below. 
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Let us estimate Q m asymptotically, as a — > 0. 
Integrate equation ([22]) over S m to get 

Qm = 2 J [f e (s),N s }ds - 2 j [N„ Bu m ]ds. (31) 

We will show in the Appendix that the second term in the right-hand side 
of the above equation is equal to — Q m plus terms negligible compared with 
\Q m \ as a — > 0. Thus, 

Q m = I [f e (s),N s ]ds, a^O. (32) 

Let us estimate the integral in the right-hand side of (|32p . 
It follows from equation f[15j) that 

[N s , f e ] = [N s , E e ] - -^[N s , [V x E e , N s }]. (33) 



If E e tends to a finite limit as a — > 0, then formula (|33|) implies that 

[N s J e ] = 0(( m ) = o(±Y a^O. (34) 

By Lemma O the operator (J — A m )~ l is bounded, so a m = O (^r), and 

Q m = O (a 2 " K ) , a -»• 0, (35) 

because the integration over S m adds factor 0(a 2 ). It will follow from our 
arguments that Q m does not vanish at almost all points, see formulas fj37[) - 

435). 

As a — )> 0, the sum (f30|) converges to the integral 

E(x) = Eq(x) + V x I g(x, y)N(y)Q(y)dy, (36) 
Jn 

where Q(y) is the function uniquely defined by the formula 

Qm = Q(x m )a 2 - K . (37) 

The function Q{y) is defined uniquely, because, as a — > the set of points 
{ x m}m=i becomes dense in The physical meaning of vector E{x) in 
equation (|36p is clear: this vector is the limit of the effective field E e (x) as 
a — > 0, and N(x) is the function from equation (I23p . 
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The function Q(y) can be expressed in terms of E: 

Q(y) = -^ L h(y)(VxE)(y). (38) 
Sou no 

This important formula is derived in Appendix. 

The factor ^ appears if D m are balls. Otherwise a tensorial factor c m , 
depending on the shape of S m , should be used in place of 

Prom equations ([5U|) and (|3"8"j) one obtains 

E{x) = E (x) - -^-V x / g(x,y)h(y)N(y)V x £7(y)dy. (39) 

Xe£ us derive physical conclusions from equation (|39|) . 
Applying the operator V x Vx to ([39]) yields 

2 , 



VxVxB = AT£ (aO 



V x (grad div - V 2 ) / g(x,y)h(y)N(y)V x E(y)dy 



n 

87T2 /* 

fc 2 S - A^tt^V x / g{x,y)h(y)N(y)V x £(y)dy 



3cL>/io 

- fc 5 

—^h(x)N(x)V x V x S(z) 
3w/io 

fc 2 £(x) - — ^-^(x)JV(a;)V x V x E 
3w/xo 

^-[V(/i(x)iV(x)),V x E{x)]. 

(40) 



Here we have used formula V x grad = 0, equation (|T0j) . and assumed for 
simplicity that h(x) is a scalar function. It follows from (|40p that 
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VxVxE = i^ 2 (x)£ - — 8 J7° , [V(%)^)),V x £(*)], (41) 



where 



^ 2 (x) = — * k 2 = cj 2 e ^. (42) 

If 

V x E = iojfi(x)H, V x H = —iue{x)E, 
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then 

VxVxB^^ijEl^VxB]. (43) 

Comparing this equation with (|4ip , one can identify the last term in (|41|) as 
coming from a variable permeability fi(x). This n(x) appears in the limiting 
medium due to the boundary currents on the surfaces S m , 1 < m < M. 
These currents appear because of the impedance boundary conditions ([7|). 

Let us identify the permeability fJ,(x). 

Denote 

:= H —h(x)N(x). 

Let e(x) = eo, eo = const, and define 

Mo 



Then K 2 = u} 2 e fi(x), and 



VfJ,(x) _ V^(x) 



Consequently, formula (|4ip has a clear physical meaning: the electromag- 
netic properties of the limiting medium are described by the variable per- 
meability: 

u ( x ) = Jf2_ = ^2 (u) 



3 Conclusions 

T/ie limiting medium is described by the new refraction coefficient K 2 (x) 
(see (|42|) ) and the new term in the equation (|4ip . This term is due to the 
spatially inhomogeneous permeability /u(x) = generated in the limiting 
medium by the boundary impedances. The field E(x) in the limiting medium 
( and in equation (|4ip ) solves equation (|39p . 

Therefore, we predict theoretically the new physical phenomenon: by em- 
bedding many small particles with suitable boundary impedances into a given 
homogeneous medium, one can create a medium with a desired spatially in- 
homogeneous permeability (|44p . 

One can create material with a desired permeability fi{x) by embedding 
small particles with suitably chosen boundary impedances. Indeed, by for- 
mula (I44p one can choose a complex-valued, in general, function h(x), and 
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a non-negative function N(x) > 0, describing the density distribution of 
the small particles, so that the right-hand side of formula (|44p will yield a 
desired function n(x). 



4 Proofs of Lemmas 

Proof of Lemma [H 

From equations (pQ) one derives (the bar stands for complex conjugate): 

/ (H -V x E - E -V xH)dx = (iufj, \H\ 2 - iue \E\ 2 )dx, 
Jd r Jd r 

where Dr := DHBr, and R > is so large that D m C Br := {x : \x\ < R} 
for all m. Recall that V • [E, H] = H-X7xE-E-X7xH. Applying the 
divergence theorem, using the radiation condition on the sphere Sr = OBr, 
and taking real part, one gets 



M ,. M 



= jr Re / [E, H] ■ Nds = Re / Cm * E t ■ E'ds, 

m=1 J Sm m=1 JSm 

where E^ is the limiting value of E l on S m from D, E l = £ m [H,N]. This 
relation and assumption (|4|) imply E^ = on S m for all m. Thus, E = H = 
in D. 

Lemma Q] is proved. □ 
Proof of Lemma [H 
If a m = A m a m , then the functions 

H = V X E , E(x)=Vx [ g(x,t)a(t)dt 

l^HO JSm 

solve equation ([1]) in D, E and H satisfy the radiation condition, and , 
condition §2$. Thus, E = H = in D. 
Consequently, 

= VxVx / g(x,t)a m (t)dt = (grad div - V 2 ) / g(x,t)a m (t)dt 
= k 2 g(x,t)a m (t)dt, x G D. 

J Sm 

This implies cr m (s) = 0. 

Lemma [2] is proved. □ 
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5 Appendix 

In Section 5.1 equation ([39]) is derived. In Section 5.2 a linear algebraic 
system (LAS) is derived for finding vectors Q m in equation ([37|) . If one 
substitutes formula ([38]) into ([36]) . then one obtains equation ([39]) . 

5.1. Derivation of the basic equation ([39]) 

Boundary condition ([7]) yields 

= [N[E e ,N]\ - ^[V x E e ,N] + [N, [V x / g(s, t)a m (t)dt, N}\ 
IUJ W Js m 

_J^[ V xVx / g(x,s)a m (t)dt,N}. 
lUfJ, Js m 

Let us denote 

f e := [N, [E e ,N]] -^[Vx i?e,iV]. 
One has V x Vx = curlcurl = graddiv — A, and 



VW g(x,t)a m (t)dt = - (V t g(x,t),a m (t)) dt = g(x, t)V t -a m (t)dt = 0, 
and 

-V^ / g(x,t)a m (t)dt = k 2 g(x,t)a m (t)dt, 

J Sm J Sm 

because —V 2 g(x,t) = k 2 g(x,t), x ^ t, see ([TO]) . Thus, using ([T3]) . one gets: 



= f e + [[ [N„ [V s g(s,t),a m (t)]]dt,N s ] + ^[a m {s),N s ] 

+ ^[N S , f g(s,t)a m (t)dt}. 

Cross multiply this by N s from the left and use the relation N s ■ a m {s) = 0, 
to obtain 

= [JV„ f e ] + [N s , [ I [N s , [W s g(s, t), u m {t)]]dt, N s ]] + ^a m (s) 
- Cm^e [iV s , [N s , / g(s,t)a m (t)dt}}. 

J Sm 



14 



Note that 

[N at [f [N s ,[X7 s g(s,t),a m (t)]]dt,N s }]= f [N s ,[W s g(s,t),a m (t)]]dt 

■J Sm J Sm 



N S ,N S ] / (N s ,[V s g(s,t),a m {t))dt), 

J Sm 

[N s ,[V s g(s,t),a m (t)])dt. 
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Consequently, 

<r m (t) = 2[f e (s),N s ) + 2( m iuje [N s , [N„ f g(s,t)a m (t)dt]] 
-2 [N s ,[V s g(s,t),a m {t))]dt := Aa m + f m , 

J Sm 



which is equation (fl4"|) . and f m := 2[f e (s), N s ], which is equation (fT5|) . 
Denote 

Qm= I (Tm(s)ds. 

J Sm 

One has 

/ [[N s , [E e (s),N s ]},N s }ds = I [E e (s), N s }ds = - f V x x E e dx. 

The term f D V x x E e dx = 0(a 3 ) is negligible compared with the terms of 
order 0(a 2 ). Let us estimate the terms of the order 0(a 2 ). One has 

/ [[V xE e ,N s ),N s ]ds = -( j VxE e ds- f iV s (V x E e ,N s )ds) 

J Sm \" Sm " Sm / 

= - / V x E e ds + -— V x E e (x m ) 

° *~>m 

8ira 2 _ „ , . 
= —VxE e (x m ), 

Here we have used the formulas 

/ V x E e ds ~ 47ra 2 V x E e (x m ), a -» 0, 

" Sm 

and 

47T(X 2 

NiNjds = —^—Sij, 
6 

15 



where S is a sphere of radius a, {Ni}f =1 are Cartesian components of the 
outer unit normal to the sphere S, and 8ij = if i ^ j, Sa = 1. 
Thus, 

0.5/ / ra (#---^( m a 2 Vx£ e M = 0(a 2 - K ), a^O (45) 
provided that 

Cm = %l, <K<1. 

a K 

Let us now show that the term J Sm Aa m ds contributes the term —Q m , so 
Q m = 0.5 I f m (s)ds(l + o(l)), a^O. (46) 

One has 



,(s,t)(N s ,a m (t)) -<j m {t)^^ ) d I 



Since 



-2 j ds j [N s ,[V s g(s,t),a m {t)]]dt 

J Sm J Srn 

= -2 / ds / di ( V s g 
= -2 I ds j dtV s g(s,t)(N s ,a m (t)) + j o m (t)dt2 I ds^^- 

J Sm Sm Sm " Sm & 

2 1 ds ^M^ = ~ 2 l rf ^,t)-i, 

J Sm s J D m 



one gets 
I : 

Therefore 



j dta m {t)2 j ds d9 ^ = - j a m {t)dt-2k 2 j dtu m {t) j dxg{a 

J Sm J Sm s J S m Sm 



I '■= —Qm + h, 
where the term I\ is negligible compared with Q rn , because 



/ dxg(x,t) = 0(a 2 ), a -> 0. 

JDm. 



Consequently, I\ is negligible compared with /. 

If f s \a m (t)\dt < oo and Q m = J s a m (t)dt ^ 0, then 



/ a m (t)dt\ > | / dta m (t) / dxg(x, 

J Sm J Sm 



t)\, 



16 



because \ J D dxg(x,t)\ = 0(a 2 ) if x G D m . If ka <C 1, then the effective 
field E e and Vx£ c change negligibly at the distances of order a, and Q m is 
proportional to a 2_K V x E e {x m ) on the surface 5 m , and therefore Q m ^ 
at all points at which V x E e {x m ) does not vanish, i.e., at almost all points 
in fi. 

One has 

| -2 J ds I dtV s g(s,t)(N s ,a m {t))\ < | j a m {t)dt\ = \Q m \, 

because \(N s ,a m (t))\ = 0(\s — t\) as \s — t\ — > 0. 
Therefore, 

f* 8* 71 / 

Qm = 0.5/ / TO (t)dt ~ Cm« 2 V X £ e M, « 0. (47) 

This yields the following formula (cf (|30p ): 
E e (x) = Eq(x) - — ^— V C m a 2 [V5(x,x m ), V x # e (aj m )], a -> 0, (48) 

m=l 

which can be written as 

£ e (x) = E (x) - -^a 2 - K V h(x m ) [V x g(x, x m ), V x E e (x m )} , (49) 

and when \x — Xj\ ~ a then the term with m = j in the sum (I49D should be 
dropped according to the definition of the effective field E e (x). 

Passing to the limit a — > in equation ([49ft . and denoting by E{x) the 
limit of the field E e {x) as a —> 0, one obtains 

£ = Eq(x) - / [S/ X g(x, y), V x E e {y)]h(y)N{y)dy, (50) 

where /i(x) is the function in the formula £ m = > an d A^(x) is the 

function in the definition of 7V(A). 

The above passage to the limit is explained in Appendix, see also [TO] . 

p. 206. This passage uses the convergence of the collocation method for 

solving equation (jM}, see [9]. Equation ([50]) for the limiting field E(x) is 

equivalent to equation ([39]) . 

5.2. Derivation of a linear algebraic system for V m 

In this Section a numerical method is developed for solving many-body 

wave scattering problem when the scatterers are small in comparison with 
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the wavelength. The method consists of a derivation of a linear algebraic 
system for finding vectors 

V m := (V x E){x m ), l<m<M. 

If V m are found, then by formulas ([55]) and (f3"T|) one finds 

Qm = -^^a 2 - K h{x m )V m , (51) 

and, by formula (f30|h the field -E(x). 

Let us derive linear algebraic system for finding V m . 

Apply Vx to equation ([15]) . let x = Xj, 1 < j < M, and replace Ylm=i 

by the sum Em^,m=i- 
Then one obtains 

8tti M 

' P 3 = 'Poj-7 ) a 2 ~ K V" {k 2 g(x,x m )h(x m )V rn +h(x m )(V m , V)V g(x,x m )}\ 2 

(52) 

where 1 < j < M, and 

P j := (V x Eo)(xj), Vj = (V x E)(xj), I < j < M. (53) 



Equation (I52p is a linear algebraic system for finding V m . In the above 
derivation we have used the known formula 

V x [A, B] = (B, V)A - (A, V)B + A(V, B) - B(V A), 

which yields 

V x [Vg, Q m ] = (Q m , V)V<? - Q m V 2 g, 

where 

3 3 d 2 g 

{Qmi V) = ^ ^ ej ^ ^ Qmi' i Qmi = {Qmj 

1=1 i'=l ' 

and ej, z = 1,2, 3, is the orthonormal Cartesian basis of M 3 . 
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